Evaluating and improving the cluster variation method entropy functional for Ising 

alloys 
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The success of the "Cluster Variation Method" (CVM) in reproducing quite accurately the free 
energies of Monte Carlo (MC) calculations on Ising models is explained in terms of identifying a 
cancellation of errors: We show that the CVM produces correlation functions that are too close 
to zero, which leads to an over-estimation of the exact energy, E, and at the same time, to an 
underestimation of —TS, so the free energy F — E — TS is more accurate than either of its parts. 
This insight explains a problem with "hybrid methods" using MC correlation functions in the CVM 
entropy expression: They give exact energies E and do not give significantly improved —TS relative 
to CVM, so they do not benefit from the above noted cancellation of errors. Additionally, "hybrid 
methods" suffer from the difficulty of adequately accounting for both ordered and disordered phases 
in a consistent way. A different technique, the "Entropic Monte Carlo" (EMC), is shown here 
to provide a means for critically evaluating the CVM entropy. Inspired by EMC results, we find 
a universal and simple correction to the CVM entropy which produces individual components of 
the free energy with MC accuracy, but is computationally much less expensive than either MC 
thermodynamic integration or EMC. 
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PACS numbers: 64.60.C, 05.50.+q 



I. INTRODUCTION 

The physics of phase transitions and phase stability of 
alloys is often couched in terms of statistical mechanics 
models on the generalized (long-range pair and multi- 
body interactions) Ising lattice and computed most ac- 
curately with Monte Carlo (MC) methodsEl. These are 
time consuming calculations and usually are thus re- 
stricted to coarse grids of chemical potential and tem- 
peratures. Further, MC simulations do not give directly 
the values of important thermodynamic variables such 
as entropy and free energy, since these quantities can- 
not be written in terms of ensemble averages. Instead, 
these are obtained laboriously by integration of thermo- 
dynamic relations from a known starting point. To rem- 
edy this situation, one often uses the less accurate Molec- 
ular Field methods, most notably the Cluster Variation 
Method (CVM) of KikuchiEl. Despite its great simplicity, 
the CVM reproduces many of the features of phase di- 
agrams obtained by the many-orders-of-magnitude more 
computer expensive MC method: For the fee nearest- 
neighbor antiferromagnetic Ising Hamiltonian with cou- 
pling constant J and zero chemical field h=0, the_transi- 
tion temperature TJ3 from MC simulations is 1.74u while 
CVM in the Tetrahedron (Tetrahedron-Octahedron) ap- 
proximation rffives 1.89 (1.81). This and other successes 
of the CVMa are even more surprising in light of the 
finding that the CVM correlation functions (the ther- 
mal average of products of the Ising spin variables) dif- 
fer considerably from the exact MC values: For exam- 
ple, in the nearest-neighbor fee antiferroniagnetic Ising 
model, the MC pair correlation functionsB at T = 1.9 



and h = are -0.208, 0.254, 0.036, 0.076 for the first 
to fourth neighbors, respectively, while the tetrahedron- 
CVM first neighbor correlation function is —0.188 and 
the tetrahedron-octahedron-CVM first and second neigh- 
bor correlation functions are —0.198 and +0.198. Thus, 
CVM correlation functions are substantially closer to 
zero (i.e., more "random") than the exact values. The 
error in tetrahedron-CVM first neighbor correlation func- 
tion leads to a ~10% error in both energy and entropy 
relative to MC (see below). However, despite such sys- 
tematic discrepancies (of ~10% or less) in reproducing 
correlation functions, the CVM seems to describe well 
thermodynamic properties (e.g., free energies) which de- 
pend on these very correlation functions. The subject of 
this paper is precisely these types of errors in CVM en- 
ergy, entropy, and free energy relative to MC. We make 
four points: 

(i) We show that a reason for the success of the CVM 
in describing the free energy is an interesting cancellation 
of errors: The closer-to-zero CVM correlation functions 
imply greater randomness and hence an over estimation of 
the internal energy compared to MC. However, the more 
random CVM correlations also lead to a larger entropy, 
and hence to an underestimation of the — TS term. Thus, 
the error in internal energy is of opposite sign to the error 
in the — TS term, so these two errors partially cancel 
in the free energy. This cancellation of errors is due to 
the fact that the CVM free energy expression may be 
obtained from a variational argument ,u but not E or S 
individually. 

(ii) Our analysis gives insight into the successes and 
failures of various approaches that attempt to improve 
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CVM by "borrowing" certain quantities from MC. In- 
deed, for some applications, one may require accuracies 
and flexibilities beyond those provided by the CVM, and 
so there is a desire in the field to find accurate "hybrid" 
methods combining the simplicity of the CVM with the 
accuracy and flexibility of MC. A natural possibility is to 
usea the correlation functions IImc (or cluster probabili- 
ties) obtained from MC simulations in the CVM expres- 
sion for free energy FqvmO^-mc) m the hope of obtain- 
ing a more accurate free energy. We demonstrate that 
these mcthodso are unlikely to succeed, as these approx- 
imations do not benefit from the cancellations of errors 
noted above. 

(iii) We use p the "Entropic Monte Carlo" (EMC) 
method of Leejj which provides a method for deter- 
mining the entropy as a function of any state variable. 
We apply EMC to the case of the CVM entropy as 
the state variable, and demonstrate that Semc('S'cvm) 
provides a means for critically evaluating the errors in 
CVM entropy. The calculation of S'emc[<S'cvm(IImc)] 
further shows that this functional accurately describes 
5 , Mc(nMc) and thus the Monte Carlo free energy. How- 
ever, this approach is computationally intensive. Finally, 
inspired by the EMC philosophy, 

(iv) We develop a functional S[Scvm(IImc)] that re- 
produces the exact SmcKHvic) very closely, and is com- 
putationally much less expensive than either MC ther- 
modynamic integration or EMC. This functional per- 
mits one to borrow from MC calculations the correlation 
functions (or cluster probabilities), evaluate the ensuing 
CVM entropy /Scvm(IImc)i an d thus obtain the nearly 
exact entropy S[Scvm(IImc)] ~ Smc(IImc) and energy 
E(Hmc)- This approach thus combines the accuracy of 
MC with the computational simplicity of the CVM. 



II. METHODS 

All of the calculations described in this paper will be 
tests of the various methods on the fee nearest-neighbor 
antiferromagnetic Ising model. 



A. CVM Quantities 

We first briefly review our notation. Let a mean a 
configuration ( "microstate" ) of Ising spins (±1) on a 
lattice. Consider a cluster ("figure") / with kf lattice 
points. The spin variable, which takes on the value 
S{(a) = — 1(+1) if there is an A (B) atom at site i of 
the figure, depends on the configuration a of spins in the 
lattice. Consider now all the clusters Rf obtained from 
the cluster / by the symmetry operations R of the space 
group of the lattice. In the CVM, we define the correla- 
tion function fl f (a) for the cluster / in the configuration 
a as the product of spin variables over the sites of /, av- 
eraged over all the figures obtained from / by the space 
group operations R: 
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The CVM treats the correlation functions {II/} as ther- 
modynamic variables. For a cluster / of kf sites there 
are 2 kf arrangements of spins ±1 at its sites. Each ar- 
rangement j has a cluster probability which is linearly 
dependent on the correlation function values for all sub- 
clusters of /. The correlation functions {11/ } (or equiva- 
lently, cluster probabilities) are determined by minimiz- 
ing the free energy, composed of the CVM internal energy 



fCF 



and the CVM entropy 



s CV M = -kY, Bf Y.p f 3 ln (pj 



(2) 



(3) 



both written as a sum over all the subclusters of the max- 
imum F. In the CVM entropy expression, one also sums 
over the arrangements of spins at the sites of each sub- 
cluster. The Barker coefficients B* can be obtained from 
purely group theoretical arguments.ErO Unless otherwise 
noted, all CVM calculations described in this paper are 
for the fee tetrahedron approximation. 

In order to examine the errors involved in the CVM, 
we first compute the accurate energy, entropy, and free 
energy from Monte Carlo simulations of the nearest- 
neighbor anti-ferromagnetic Ising model at h=0. 



B. Monte Carlo Quantities 

A Monte Carlo cell of 1728 sites was used with 10 6 
Monte Carlo steps per site at each temperature. Al- 
though finite-size effects were not taken into account, 
the calculated heat capacity showed a sharp peak at 
T = 1.77, within ~l-2% of the most precise values for 
the transition temperature given in the literature ~1.74- 
1.75.H The energy is given directly from MC, while the 
entropy is obtained by thermodynamic integration down 
from infinite temperature: 



S(T) = k\n2- 



E(T) 
T 



E(T)d(l/T) (4) 



(The entropy was also obtained by integrating the heat 
capacity down in temperature, however, this method was 
found to be less efficient in that it required a finer grid 
of temperatures near the transition for equal accuracy). 

The correlation functions IImc for a figure / were ob- 
tained by taking the thermal average (over the 10 6 Monte 
Carlo steps) of the product of Ising spin variables over 
the sites 1,2, ...,kf of all symmetry-equivalent figures / 
[i.e., the thermal average of Eq. (|l|)]. 
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FIG. 1. Energy, entropy, and free energy as a function of 
temperature for the nearest-neighbor anti-ferromagnetic Ising 
model. All quantities are given in dimensionless units: ksT/J 
for temperature, energies are given normalized by J, and en- 
tropies are given normalized by fc_g. (a) Results obtained from 
Monte Carlo simulations and thermodynamic integration, (b) 
Errors in standard CVM (disordered phase symmetry) com- 
pared to Monte Carlo. We have used only the CVM entropy 
expression with the disordered phase symmetry. Thus, in 
(b), differences with Monte Carlo for temperatures below the 
CVM transition (T=1.89) are overestimated and hence are 
shown as dashed lines (see text). 



III. RESULTS 

A. Analysis of CVM errors vis-a-vis MC simulations 

The energy, entropy, and free energy obtained from 
Monte Carlo simulations are shown in Fig. [j]a. The first- 
order transition at T ~ 1.77 is evident from the discon- 
tinuity in energy and entropy. We have also computed 
the energy, entropy, and free energy predicted by CVM 
(in the tetrahedron approximation) . By comparing these 
CVM results with the "exact" Monte Carlo results in Fig. 
|l|a, we may ascertain the errors in thermodynamic func- 
tions of the CVM. The differences 5E, S(-TS), and SF 
between the respective CVM and Monte Carlo functions 
are shown in Fig. [l]b. In our CVM calculations, we have 
used only the CVM entropy expression with the disor- 
dered phase symmetry. Thus, in Fig. |l|b differences with 



Monte Carlo for temperatures below the CVM transition 
(T=1.89) arc shown as dashed lines. Obviously, practi- 
tioners of CVM would correctly impose a lower symmetry 
on the entropy expression below the transition to the or- 
dered phase and would not use the disordered phase sym- 
metry in this temperature range. The reason we use the 
CVM disordered phase symmetry down to low tempera- 
ture is due to our wish to combine CVM methods with 
MC, which as we describe in the next section, is problem- 
atic when using ordered expressions for CVM entropy. 

As indicated in the Introduction, the CVM correlation 
functions {II/} are closer to zero than the MC values. By 
Eq. (||) , the CVM internal energy is less negative relative 
to Monte Carlo (thus, SE > in Fig. [l]b), demonstrat- 
ing that the energetic effect of short-range order in CVM 
is underestimated, and hence the CVM internal energy 
is "more random" than that of Monte Carlo. The en- 
tropy of CVM is overestimated (6(—TS) < 0) relative to 
Monte Carlo, again indicating a more random solution 
than Monte Carlo. However, the error in the CVM free 
energy SF is simply the sum of the errors SE + S(—TS). 
Since 5E and S(—TS) have opposite sign, they partially 
cancel, and give an error in free energy which is consider- 
ably smaller in magnitude than either the error in energy 
or in entropy Thus, owing to the variational nature of 
the CVM,H the free energy of CVM is more accurate than 
one might expect from considering cither the energy or 
entropy alone. 



B. Using MC correlation functions in CVM 
calculations: Absence of cancellation of errors 

Our foregoing discussion sheds light on a hybrid 
method which, naively thinking, might combine the ac- 
curacy of Monte Carlo with the simplicity of CVM. In 
this method, one uses the correlation functions {II/} (or 
equivalently, the cluster probabilities {Pj}) of MC sim- 
ulations in the expressions for CVM entropy [Eq. (||)] 
and energy [Eq. (g)]. We refer to this method as the 
"-FcvmCIImc)" method. This method would, of course, 
require one to perform a Monte Carlo simulation for each 
composition and temperature of interest; however, one 
could, in principle, obtain the entropy at each point from 
a single Monte Carlo simulation (i.e., one composition 
and one temperature) rather than a series of Monte Carlo 
calculations which would be required for thermodynamic 
integration of the entropy. Since the Monte Carlo cor- 
relation functions are used in this method in Eq. (|^), 
there is no error in energy (SE = 0). Thus, the error 
in free energy, shown in Fig. |^, is equal to the error in 
entropy: SF = 5(—TS). Since there is no error in en- 
ergy in this method, there is no cancellation of errors. 
Hence, even though the exact Monte Carlo correlation 
functions are used in the "Fcvm(J1mc)" method, it pro- 
duces less accurate free energies than standard CVM: For 
example, at T=1.92, the free energies as given by Monte 
Carlo, "F C vM(n M c)", and CVM, are = -2.109, -2.050, 
and -2.067, respectively. (For comparison, CVM in the 
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FIG. 2. Error in free energy as a function of temperature 
for the nearest-neighbor anti-ferromagnetic Ising model ob- 
tained using the "FcvmOImc)'' method. All quantities are 
given in dimensionless units: ksT/J for temperature, ener- 
gies are given normalized by J, and entropies are given nor- 
malized by ks- In the "Jcvm(IImc)" method, MC correlation 
functions and cluster probabilities are used in the CVM ex- 
pressions for energy and entropy, respectively. We have used 
only the CVM entropy expression with the disordered phase 
symmetry. Thus, differences with Monte Carlo for tempera- 
tures below the CVM transition (T=1.89) are overestimated 
and hence are shown as dashed lines (see text). Note that in 
this method, the energy is precisely that of Monte Carlo, thus 
the error in — TS is also the error in free energy. 



Tetrahedron-Octahedron approximation gives F — -2.094 
for this temperature.) 

Other disadvantages of the -Fcvm(Hmc) approach are 
illustrated in Fig. 0, showing a comparison of the en- 
tropies as a function of temperature as calculated by 
MC (Sue), by standard CVM [ScvmCIIcvm)], and by 
the Fcvm(IImc) method [Scvm(IImc)]- (The othe r two 
curves of Fig. are discussed in Sections III C and [II D 
below.) One can again see that (i) standard CVM (solid 
line) overestimates the entropy at high temperatures rel- 
ative to Monte Carlo (open squares), (ii) the CVM en- 
tropy of the disordered phase is not applicable at low 
temperatures, and (hi) the Fcvm(^uc) method under- 
estimates the entropy at high temperatures, and at low 
temperatures this entropy takes on the unphysical values 
ScvM(nMc) < 0. These unphysical values are a result of 
the fact that the expression for the disordered Scv m al- 
lows negative values for many atomic configurations. For 
instance, in the fee lattice, the simplest ordered configu- 
rations such as Llo, ili, LI2 have negative Scvm val- 
ues. (These negative values of CVM entropy are likely to 
persist no matter what sized maximal cluster is used, and 
thus will always lead to difficulties with the -FcvmCIImc) 
method at lower temperatures, since the CVM entropy 
will incorrectly tend to negative values rather than zero.) 
Of course, one might argue that these ordered configu- 
rations only possess negative CVM entropy when evalu- 
ated with the CVM expression for the disordered phase, 
whereas when evaluated with the CVM expressions for 
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FIG. 3. Entropy versus temperature for the near- 
est-neighbor anti-ferromagnetic Ising model. The open 
squares (connected by a solid line) is the result of Monte Carlo 
simulations, the solid line is standard CVM, the long dashed 
line is the CVM entropy expression evaluated with the Monte 
Carlo cluster probabilities, the thin short dashed line is the 
result of the entropic Monte Carlo calculations, and the thick 
short dashed line is the simple correction to the CVM, the 
"modified CVM". 



the corresponding ordered phases, they will have non- 
negative entropies. However, this illuminates another po- 
tential problem with the FcvmO^mc) method: In Monte 
Carlo simulations, the presence of anti-phase boundaries 
and finite-sized domains of long-range order preclude 
one from unambiguously defining the distinction between 
sublattices and extent of long-range order present in the 
simulation. However, an ordered CVM expression for the 
entropy is written in terms of correlation functions and 
cluster probabilities for all the symmetry-distinct figures 
for the symmetry of the long-range ordered phase. Thus, 
the CVM ordered entropy expressions presuppose the do- 
mains of long-range order are infinite in size, and hence 
the distinction of various sublattices in the ordered phase 
is unambiguous. Thus, using the "i^cvM^Mc)" method 
with ordered CVM entropy expressions is not practical 
because one doesn't know from the MC simulations pre- 
cisely how to divide the MC simulations into sublattices 
of long-range order and hence one does not even know 
from the MC simulations which ordered CVM expres- 
sion to use. Also one does not know at what tempera- 
ture to change the symmetry of the CVM to the ordered 
entropy expression. Thus, an ideal method combining 
Monte Carlo and CVM would only use a single expression 
(e.g., the disordered CVM expression) for the entropy at 
both low and high temperatures. We next describe such 
a method called "Entropic Monte Carlo" (EMC). 



C. The entropic Monte Carlo method: A critical 
evaluation of CVM entropy 

LeeB has shown a practical way to determine the en- 
tropy of a Monte Carlo cell as a function of any state vari- 
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able. We call this method Entropic Monte Carlo (EMC). 
Though in his paper Lee applies the EMC method to the 
case of the energy as state variable in a quantized system, 
here we describe instead 5emc in terms of the state vari- 
able Scvm(&) which is a continuous, not quantized, vari- 
able. Our strategy will be to calculate Semc (Scvm) by 
the method of Lee, and then insert Scvm(JImc) into this 
expression, giving Semc[Scvm(^mc)} which we write as 
-Semc(IImc)- We will show that this function reproduces 
very well Smc- First, we describe how S-emc(Scvm) is 
calculated: 

The EMC method is a self-consistent process in which 
each iteration is made from a series of Monte Carlo 
sweeps where the driving "energy" E(a) of the Monte 
Carlo equations is not the true energy contained in the 
sample but an approximation to the entropy: 



E(a) = E[S C vm{(t)\ 



(5) 



which depends on the configuration a through the func- 
tion E[Scvm], whose argument is the CVM entropy (per 
site) calculated with Eq.(|^) for the cluster probabilities 



p. (a) of the configuration a. The function E [Scvm] is 
assumed to be monotonic. The EMC dynamics are given 
by the detailed balance condition 



exp [-E (a,)} W(i - j) = cxp [-E (a,)} W(j 



(6) 



where W are transition rates and E(a) is given by Eq. 
(|J). After many MC sweeps of the full lattice, one obtains 
a histogram 



HIS, 



CVM 



Xd (Scvm) exp [-E (Scvm)] (7) 



where d {Scvm) is the number of configurations with a 
given value Scvm of the CVM entropy (degeneracy in 
Scvm) an d X is a constant of proportionality. We dis- 
tinguish Scvm, which is a numerical argument attaining 
certain value, from Scvm which is a function of both 
the microstate [through Eq.(||)] and the cluster probabil- 
ities pj (a) . The density of states function of the CVM 
entropy Scvm is given by 

D(Scvm) — ^ d (Scvm) 8 (Scvm — Scvm (c)) 

= (1/V) exp (E (Scvm)) H (S C vm) x 

S (Scvm — Scvm (&)) (8) 
and the entropy Semc — $(§cvm) per site is denned as 

/Scvm 
D(£)d£ (9) 
-oo 

where V is the number of sites in the MC cell. 

From Eqs.(§) and Eq.(|), we may obtain the difference 
between the entropy at two different values Scvm- 

S{Scvm) ~ S(Scvm) = 



eMNE[S C vM(<y)}) x 



f, Scvm(t)<S { ^ } vm 



H[S C vm((t)]} 

4 m £ 



exp(NE[S C vM(o-)}) x 



<?, Scvm(<?)<ScI 



H[S C vm(<j)}} 



(10) 



which is the basic equation used to determine the entropy 
from the EMC runs. On the right-hand side of Eq. (px|), 
the sums are over the microstates a obtained in the EMC 
sweeps whose CVM entropy Scvm(o~) are smaller than 

Scvm or ^cvm- As pointed out by Lee0, the entropy 
determination becomes especially simple when the inter- 
action E(Scvm) is such that the histogram H [Scvm] 
is uniformly distributed and has little dependence on 
Scvm- In this case, Eq. (|10|) becomes 



S{S C y M ) 



S(S CVM ) 



N 



ln[ 



E 



a, Scvm(c)<S i ^ ) vm 



E 



eMNE[S CV M(o-)])] 
exp(NE[ScvM(<r)]'. 



(11) 



Scvm(<t)<S^ c 



Because of the factor of N, the exponential in Eq. (|1]^ 
is a rapidly increasing function of Scvm, and hence only 
the extremes contribute significantly to the sums, or 

S(Scvm) ~ S(S CVM ) ~ JtE(Scvm) ~ JjE(Scvm) 

(12) 

This equation then suggests the self-consistent procedure 
for determining the entropy: From a crude estimate of 
S (Scvm), we use Eq. (|12|) to obtain the interaction 
E (Scvm) with which we make EMC runs, with which we 
recalculate S (Scvm) from the basic Eq. (|To|) . This pro- 
cess is taken to self-consistency. When self-consistency 
is reached, (i) the histogram H[Scvm] (the number of 
microstates in each small range of Scvm obtained in a 
series of MC sweeps) is nearly constant, independent of 
the value of Scvm and (ii) the driving "energy" which 
is the approximation to the entropy [Eq. (j|)], becomes 
equal to the entropy calculated from the density of states 
[Eq. (g)], or in other words, the EMC entropy becomes 
exact. An important aspect of the EMC is that the cal- 
culated entropy functional form of S(Scvm) does not 
depend on any particular Ising Hamiltonian, (so long as 
the important correlations are contained within the CVM 
maximal cluster), because the role of the "energy" driv- 
ing the EMC calculations is played in Eq. (|o|) by the 
entropy itself. 

Figure |] shows a typical result Semc(Scvm) of the 
EMC calculations, using a MC cell with N = 12 3 = 
1728 sites. We also performed EMC calculations with 
different MC cells, with the results being slightly different 
for the larger negative values of Scvm- The negative 
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FIG. 4. Entropic 

Monte Carlo results for Semc = S(Scvm)- Filled circles are 
the EMC calculations, and the solid line is the EMC-inspired 
entropy functional S(Scvm)- The dashed line is line of unit 
slope simply to guide the eye. EMC was performed for a cell 
of 12 3 = 1728 sites, using the CVM tetrahedron expression for 
the disordered entropy. Note that many configurations cor- 
respond to negative CVM entropy, with the most negative 
(for all configurations with <16 atoms per cell) being the Llrj 
configuration. 



values of Scvm correspond to configurations of atoms 
with higher symmetry, usually associated with smaller 
repeat units, thus explaining why the curve depends to 
some extent on the size and shape of the MC cell for 
the negative values of Scvm- For instance, for an EMC 
cell with an odd number of sites (e.g., II 3 ), one could 
never obtain the stoichiometric configuration LIq with 
its large negative CVM entropy. In fact, simple high- 
symmetry configurations such as £lrj, Ll±, LI2 all have 
negative values of Scvm- Examining the CVM entcppy 
for all configurations with up to 16 atoms per cell,E£l we 
found the configurations with the most negative CVM 
entropy had very small unit cells. The largest negative 
CVM entropy occurs for LIq for which Scvm (LIq) = 
-1.34 ln2. 

The results Semc[Scvm(Hmc)] 01 EMC are shown in 
Fig. ||, where they are contrasted with—the results of 
Monte Carlo, CVM, and "FcvM(n M c)'' B By compar- 
ing S E Mc(ilMc) and S C vm(H C vm) with Sue, we see that 
the EMC and CVM entropies are equally accurate at high 
temperatures. Remarkably however, the EMC method 
also produces extremely accurate entropies at low tem- 
peratures, in qualitative contrast with the "-Fcvm(Hmc)" 
method. Thus, even though one only uses a single disor- 
dered expression for the CVM entropy in the EMC cal- 
culations, the EMC reproduces both high temperature 
(disordered) and low temperature (ordered) entropy val- 
ues, with no need to change the CVM entropy expression 
at any point. Although the internal energy in EMC is ex- 
act (so this method does not b enefit from the cancellation 
of errors noted in Sec. [II A for the CVM), we see that 
EMC does not need to be correct due to cancellation of 
errors. Instead, it is accurate because its individual terms 



(E and —TS) are accurate. 

The EMC, like standard Monte Carlo, can be a compu- 
tationally laborious procedure. However, our EMC cal- 
culations of S (Scvm) suggest a very simple functional 
S(Scvm) which is appealing because the correction does 
not require one to perform an EMC calculation. We next 
describe this simple correction. 



D. An EMC-inspired new entropy functional 

While ScvM(nMc) 
can be inaccurate, Semc [Scvm (IImc)] is accurate but 
computationally expensive. Thus we will now develop 
a new functional S[Scvm(Hmc)] which is both accurate 
and inexpensive. 

The EMC results of Fig. [| permit one to guess the 
behavior of the "exact" entropy S (Scvm) (hi the limit of 
N — > 00) as a function of the CVM entropy obtained from 
a "good" maximal cluster (e.g., the tetrahedron or the 
tetrahedron-octahedron). This "true" entropy function 
S(Scvm) should have the following properties: 

(i) The most positive Scvm (2) entropy S^y^(x) 



should correspond to the exact entropy for this case, 
i.e., the ideal mixing entropy: S(S^y^) = S°(x) = 
— fcs[cclnx + (1 — x) ln(l — x)]. 

(ii) The slope of S(Scvm) at the maximum value of 
S'cVM = S° should be unity because for nearly ran- 
dom configurations the CVM approaches the exact re- 
sult: ,„ dS I go = 1. 
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(iii) The most negative value of the CVM entropy 
should correspond to zero "true" entropy. Thus, 
S(Scvm) = ^' ^ nc conn g ur ation with most negative 
CVM entropy can be found by examining all configura- 
tions up to some maximum unit-cell size, as described in 
Rcf. [lO. For instance, for the tetrahedron CVM, Ll has 
the most negative CVM entropy [Scvm = Scvm (Llo) = 
— 1.341n2]. This point is indicated in Fig. ||by a square. 

(iv) The function S (Scvm) should increase monoton- 
ically with Scvm as one can see from Eq. (|l0|). Also, 
S (Scvm) has a positive curvature [due to the exponent 
of N in the right-hand side of Eq. (fij)]. 

We select a simple functional form for S (Scvm) which 
satisfies (i)-(iv) above (but otherwise possesses no spe- 
cial physical meansing.) However, use of this simple 
form will provide a means of evaluating energies and en- 
tropies which (a) is computationally much more efficient 
than either MC thermodynamic integration or EMC, (b) 
possesses MC accuracy, and (c) may be extended to 
use any maximal cluster of the CVM. The functional 
form we choose for the approximation S(Scvm) to the 
true S (Scvm) which satisfies the four properties (i)-(iv) 
above is: 

S(Scvm) = (S° - — ) + — expM^i - 1)] (13) 
where a is the solution of 

cO cO cMIN 

= (5°-^) + ^expK%t-l)]. (14) 
a a S u 



G 



In the case of the tetrahedron CVM, a = 0.86917. The 
function in Eq. ( p~4|) depends only on a single parameter 
^CVM which can be easily estimated, from an enumera- 
tion of small- unit-cell configurationstj using any maximal 
cluster of the CVM method. (In contr ast, an EMC calcu- 
lation as we described in Section III C is only practical for 
the tetrahedron CVM approximation.) While the nega- 
tive Scvm configurations, which correspond to highly 
symmetric arrangements of atoms, have no meaning in 
the standard CVM procedure (since in standard CVM 
one would use a different expression for CVM entropy to 
describe ordered phases), the main merit of a correction 
such as Eqs. (|l3|)-([l4|) is to restore these highly ordered 
configurations into a single CVM expression by attribut- 
ing a non-negative entropy to them. Naturally this cor- 
rection, when used together with the Monte Carlo cor- 
relation functions, will be especially important near the 
transition when the ordered configurations begin to be 
important. 

To test these ideas, we calculated S[Scvm(IImc)] = 
S(n MC ), where S is given by Eq. @, and the CVM 
is executed within the tetrahedron approximation. One 
sees (Fig. ||) that this approach presents a remarkable im- 
provement over the Fcvm(IImc) method in all tempera- 
ture ranges, especially below the transition temperature 
where Scvm(IImc) is negative. Also, the simple func- 
tional represented by S(IImc) effectively retains all of 
the improvements over Fcvm(IImc) that were obtained 
by the full EMC calculation. In fact, for high temper- 
atures, the S(Scvm) approach is even closer to the ex- 
act Monte Carlo results than the EMC calculations on 
which it was based! This fact can be understood by 
examining the EMC calculations in Fig. |[ Figure |] 
shows that the EMC calculations using a cell of 1728 
sites do not reproduce property (ii) above, — [go = 1. 
In fact, the slope of the EMC curve in Fig. |y is only 
about 0.82 at the maximum value of the entropy. EMC 
simulations with even smaller cells were typically found 
to possess even smaller slopes. Presumably, for larger 
EMC simulations one would approach the correct slope 
of unity. The slope of EMC being smaller than unity 
means that as one comes down from infinite temperature, 
or Scvm = In 2, the EMC entropy maintains a larger 
value than it should. This explains why SemC > Smc 
for temperatures above the critical temperature. Because 
the S(Scvm) a PP roac h was constructed to obey the re- 



quirement 



dS 
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of EMC caused by the finite-sized simulation cell, and 
hence improves the entropy above the transition. The 
reason that the derivative of EMC is less than one for 
small cell sizes is due to the negative Scvm configura- 
tions, which are typically high-symmetry, small-unit-cell 
configurations. Thus, these negative Scvm states are 
represented more in small EMC cells relative to config- 
urations with large-unit-cells and low symmetry. If the 
density of states in Eq. (^) becomes artificially large for 
the negative region of Scvm (due to small EMC cells), 
it will have to be compensated by an artificially small 
density of states in the region of positive Scvm [since 



| go = 1, it corrects the error in the slope 



the integral in Eq. (|9|) is constrained by the fact that 
it must be 2 N at Scvm = In 2.] Thus, the integral will 
grow more slowly than it should for CVM entropies ap- 
proaching In 2, and hence the slope will be less than one. 



IV. SUMMARY 



The main accomplishment of this paper is to suggest 
a simple functional 5 (Scvm) [Eqs. (|l^)-([l^)1 that im- 
proves the CVM entropy. The development was based 
on insights gained from our analysis of the CVM free 
energy (which showed cancellation of energetig vs. en- 
tropic errors), and from the EMC philosophyQ permit- 
ting one to express the true entropy as a functional of an 
approximate, but deterministic entropy. The new func- 
tional S(Scvm) can be used in future applications either 
with CVM alone [simply by replacing the CVM entropy 
with Eq. ([p?]) in any existing CVM program], or with a 
combination of CVM and IImc (using S[Scvm(IImc)] as 
described in this paper). 
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